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Abstract 

We discuss asymptotic properties of a family of discrete probability measures 
which may be used to model particle configurations with a wall on a set of discrete 
nodes. The correlations are shown to be determinantal and are expressed in terms of 
discrete orthogonal polynomials. As an application we study random tilings of the 
half-hexagon or, equivalently, configurations of non-intersecting lattice paths above 
a wall, so called water melons. 

1 Introduction 

In this article we are concerned with discrete probability distributions which are closely 
related to discrete orthogonal polynomial ensembles. Given a positive weight function 
defined on a set of nodes X^ — {%n,o < ■ ■ ■ < %n,n-i} contained in an interval [a, b], the 
probability distribution on the set of fc-tuples (x\, . . . , Xk) E X^ with x\ < X2 < ■ ■ ■ < Xk 
(k < N), given by 

1 k 

p (N ' k \x 1 ,X 2 , ■ ■ • ,X k ) = — Y\ W N(%i) 1 f (Xj-Xi) 2 (1.1) 

1=1 l<i<j<k 

is called discrete orthogonal polynomial ensemble. Here Zk is a normalisation constant. We 
refer to this ensemble as DOPE(N, k). p <yN,k \x\ 1 . . . , Xk) can be viewed as the probability 
of finding a configuration of k particles located at the sites Xi, . . . ,Xk E X N . Continuous 
siblings of these ensembles describe the joint distribution of eigenvalues of a random matrix 
from the Gaussian unitary ensemble [5j HO], or the joint distribution of k non-intersecting 
Brownian bridges, see [6]. 



1 



The results of this paper deal with a related ensemble involving an even weight function 
w 2 n (i-e. w 2 n{ x ) = u , 2n(— x)) which is defined on a symmetric set of nodes X 2N = 
Y N U -Y N C [-b, b], where Y N = {x 2NjN , x 2N , 2N -i} Q [0, 6]. We define DOPE syra (N, k) 
to be the ensemble of fc-tuples (xi, . . . , x^) E Y N , x\ < ... < Xk, equipped with the 
probability distribution 

1 h 2 

pS£ ) (x 1 ,x 2 ,...,x k ) = Y[xfw 2N (xi) Yl ( x2 j- x i) ' ( L2 ) 

k i=i i<«<i<fc 

where ^ ym is a normalisation constant. This ensemble may serve as a model for particle 
configurations in the presence of an impenetrable wall, as e.g. in the half-hexagon tilings 
discussed below. These in turn are studied in connection with an ensemble of random 
matrices in [3]. 

In Section |2| it is shown that the ensembles DOPE sym (N,k) and DOPE(2N,2k) with 
the even weight function w 2 n are closely related via their correlation functions. In the 
monograph pQ asymptotic results on discrete orthogonal polynomial ensembles are estab- 
lished for the limit k, N — > oo with k ~ cN. These are applied to prove similar results for 
the DOPE sym (N, k) ensemble. 

Remark: Throughout the paper, when dealing with the ensemble DOPE(2N,2k) we 
assume that w 2 n is an even function and the set of nodes X 2 ^ is of the above symmetric 
form. 



2 Correlation functions 



We first summarise some facts on the DOPE(N, k) ensembles (1.1 ), the most important of 
which is the determinantal structure of the correlation functions. It involves the eponymous 
orthogonal polynomials. A scalar product on the set of complex valued functions on 
is associated to the weight function wn via 

N-l 

(f, g) ^ Wn fcw) f t x *,i) 9 (x N>i ). (2.1) 

i=0 

By applying the Gram-Schmidt procedure to the sequence of monomials 1, x, . . . i* -1 we 
obtain a family of orthonormal polynomials Pn,o, ■ ■ ■ iPn,n-u i- e - t ne degree of p^j is equal 
to j and the relation 

N-l 

} u w N {x N ^p Ntk (x N>i )p N) i (x N>i ) = S ki 

t=0 

holds. Note that since the nodes, the weights and the coefficients of the PN,k are rea l; we 
can omit complex conjugation. Furthermore the leading coefficient 7^ of pN,k is assumed 
to be positive. Denote by 7Tjv,fe := ^\pN,k the kth. monic orthogonal polynomial. The m- 
point correlation functions R^ k (xi, . . . ,x m ) describe the probability that a configuration 
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of k particles in DOPE(N, k) contains particles at each of the m sites x±, . . . , x m (m < k). 
In particular the one-point correlation function k (x) equals the probability of finding a 
particle at x. We have 



(2.2) 



Rff k (xi, . . . , x m ) = P(particles at each of the sites xi, . . . , x m ) 
= det (K N)k (xi,Xj)) i j=K m , 

where for x,y G the correlation kernel K^^ k (x,y) is given by 

fe-i 

K N>k (x,y) = y / w N (x)w N (y) y^ y PN,n{x)pN,n{v) 

n=0 

I 1—\ T\ lN,k-l PN,k( X )PN,k-l(y) - PN,k(v)PN,k-l(x) 

= \fw N {x)w N {y) 

lN,k x-y 

if x ^ y, and otherwise 

K N)k (x,x) = w N {x) ■ lN ' k ~ l . (j/ Nk ( x )p Ntk -i(x) ~ p' N , k -x( x )PN,k(x)) . (2.4) 

lN,k 



(2.3) 



The second "=" in equation (2.3) is known as Christoffel-Darboux formula. The derivations 
of both the particular form of the correlation functions and the summation formula are 
fundamental calculations in random matrix theory, see (5j Ch. 4] or [TQl Ch. 5]. 

We can also obtain determinantal representations for the m-point correlation functions 
of the DOPE sym (N,k) ensembles (1.2). To this end we need monic polynomials qj(z) of 
degree j, j = 0, . . . , N — 1 with the property 

^2 x2w 2N(x)qi (x 2 ) qj (x 2 ) = %. 

Once these are at hand, we can mimic the standard random matrix computations [5J Ch. 
4] and find a determinantal representation of the m-point correlation function with kernel 

fe-i 

k n™( x > v) = V x2w 2n(x) Vy 2 w 2N (y) e * qn qn (^) • 

n=0 

Peter Forrester pointed out how these polynomials qj can be obtained. Consider the monic 
polynomials orthogonal with respect to the even weight function w 2 n on the set of nodes 
X 2N = Y N U -Y N , 

22 w 2N( x )nw,i ( x ) K 2 N,j (x) = — 2 Sij. 

XGX 2N 

Since w 2 n is even, it follows easily from the Gram-Schmidt procedure that iT2N,2j(x) is even 
and ir 2 N,2j+i(x) is odd. In particular we have ir 2 N,2j+i(x) = xqj (x 2 ) for a monic polynomial 
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qj of degree j. The qj satisfy 



} u x 2 w N {x)qi (x 2 ) qj (x 2 ) 
xeY N 

} u W N (x)7T 2 N,2i+l ( x ) K2N,2j+l ( x ) 
x£Y N 

9 ^2 W N( x )^2N,2i+l( x )^2N,2j+l( x ) 



(2.5) 



2 7 | 



-Si 



xdX 2N ~ '2JV,2i+l 

and hence are the sought for polynomials. For x, y £ Y/v the correlation kernel for the 



ensemble (1.2) can be written as 



k-1 



K nJ( X i V) = V x2w 2N(x) \fy 2 W 2N {y) 2 72iV,2n+l<?n (x) q n (y) 

n=0 

k-1 

= \/w 2N (x) \J W 2N {y) Y 2 727V,2n+l 7r 2iV,2n+l ( x )^2N,2n+l (v) 
n=0 
k-1 

=2 a/ W 2 j\f ( x ) \/ W 2 N (y) ^ P™,2n+1 ( x )p2N ,2n+l (y) 



(2.6) 



■\/w 2N (x)^w 2N (y) 



n=0 
~2k-l 



2fc-l 



^2p2N,n( X )P2N,n(y) ~ ^ P2N,n( X )P2N,n(~y) 



n=0 



n=0 



--K 2Ny2k (x, y) - K 2Nj2k (x, -y). 



In the following we will show that for x,y £ Yjv and x, y > e > the summand 
K 2 N,2k{x, —y) tends to zero in the considered limit, and hence the correlation kernels 
K 2 N,2k and asymptotically coincide. For nodes which are in 0(1/ N) distance to 0, 

this is in general not the case. 

Remark, i) The DOPE sym (N,k) ensembles can also be defined with a weight w 2 n+i on 
a symmetric set of nodes X 2 ^ + i = — Y/v U {0}U Ypj. Configurations containing occur with 



probability and the calculations of equation (2.5) also apply in this situation, leading to 



a kernel = K 2N+12k (x,y) — K 2N+ i t2k (x, —y). For notational convenience, we discuss 
the case of an even number of nodes. 

ii) A similar relation is present between ensembles of non-intersecting Brownian motions 
and excursions [33], the continuous analogues of the non-intersecting lattice path models 
corresponding to random tilings of the full and the half-hexagon, respectively. 

3 Discrete orthogonal polynomials 

In order to state our asymptotic results, we review some terminology and results from the 
monograph [1] . In applications as the arctic circle phenomenon, the asymptotic behaviour 
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of PN,k{z) and i"Gv,& as N and k simultaneously tend to infinity plays a crucial role. The 
results of pQ are obtained under some technical assumptions [TJ Section 1.2] on the weight, 
nodes and number of particles. These are in particular fulfilled in the situation of the 
(half-) hexagon tilings. 

3.1 Basic assumptions 

The nodes 

We assume the existence of a node density function p°, which is real- analytic in a complex 
neighbourhood of [a, b], strictly positive in [a, b] and satisfies a normalisation condition and 
a certain quantisation rule, namely 

/ p°(x)dx = 1 and / p°(x)dx = — ±— , n = 0, . . . , N - 1. (3.1) 

J a J a 

The weight function 

We assume that we can write the weight function in the form 

JV-l 



W N 



{x N ,n) = (-lf-^ e - NVN M H ( XN>n - x^y 1 , (3.2) 



m=0 



where Vn(x) is a real-analytic function defined in a neighbourhood G of [a, b]. Furthermore 

V N (x) = V(x) + (3.3) 

where V(x) is a fixed real-analytic potential function independent of iV and 

limsupsup ^(z) | < 00. 

Af-s>oo zeG 

As opposed to V(x), the correction rj(x) may depend on N. 
The number of particles 

The number of particles k (at the same time the degree of the involved polynomials) and 
the number of nodes iV are related by 

k = cN + k, (3.4) 

where c G (0,1) and k remains bounded as iV — > 00. 

Further assumptions are difficult to express explicitly in terms of the nodes and the 



weight and postponed to Section 3.3 
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3.2 The associated equilibrium energy problem 

The asymptotic behaviour of n N<k is closely related to the asymptotic distribution of zeroes 
in the interval [a, b], wich in turn can be expressed in terms of quantities arising in a related 
constrained variational problem [9] . Define a real-analytic function (p by 

f b 

<p(x):=V(x)+ / log|x — y\dy. (3.5) 

J a 



Note that according to the representation (3.2) wn(x) is asymptotically equal to 

w N (x) ~ e -"K*)-»K*). (3.6) 



Recall that T^N,k minimises the scalar product (2.1) among all monic polynomials of 
degree k. This relates the asymptotic distribution of zeroes of tt^^ to (the variational 
derivative of) the quadratic functional of Borel measures on [a, b] defined by 

E c [fi,V] := E c \p] := c / log - - dfi(x)dfi(y) + / tp(x)dfi(x). (3.7) 

Ja Ja V\ Ja 

We are looking for a measure yU^ in which minimises E c \p\ subject to the normalisation 
condition 

rb 

dfx(x) = 1 (3.8) 



and the upper and lower constraints 

< j dfi(x) < - j p°(x)dx for every a Borel set B in [a, b]. (3.9) 

We refer to /^ iT1 as the equilibrium measure. The constraints are due to the fact that all 
zeroes of 7r N)k are contained in the interval [xn,o, x n,n-i] an d a closed interval [xn iTI , x^ tn+ i\ 
between two consecutive nodes contains at most one zero of 7rjv,fc- Note that minimising 



E c [jj\ simply subject to the normalisation condition (3.8) is formally like seeking a critical 
point of 



F c [ii] = E c [fi] -l c dfi(x 

J a 



with a Lagrange multiplier l c := l c [V}. When fi = /i^in' is a real constant. /x^ in is known 
to be unique and it has a piecewise analytic density dfj,'^ lin (x)/dx. Points of non-analyticity 
are finite in number and do not occur in points where the upper and lower constraints 
<i/i^ in (x) / 'dx > and d^ in (x) j ' dx < 1/c hold strictly and simultaneously. 

3.3 The equilibrium measure: related quantities and further as- 
sumptions: 

The constraints give rise to the following 
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Definition 3.1. A band is a maximal open subinterval of [a,b] where /i min is a measure 
with a real- analytic density d^ in (x)/dx which satisfies < d^ in (x) / dx < p°(x)/c. A 
void is a maximal open subinterval of [a,b] in which d^i^ oin (x)/dx = 0, i.e. meets the 
lower constraint. A saturated region is a maximal open subinterval of [a, b] in which 
dfi^^x) / dx = 1/c and hence meets the upper constraint. If no stress is put on the active 
constraint, voids and saturated regions are referred to as gaps. 



As announced in Section 3.1 we make some further assumptions on the weight and the 
nodes which are expressed in terms of the equilibrium measure, cf. [TJ Section 2.1.2]. Let 
T be the closed set of points x G [a, b], where d^ in (x) / dx = p°(x)/c or d/J^^x) / dx = 0, 



i.e. one of the constraints (3.9) is active. We assume that the connected components of J 7 
have non-empty interior and a, b G J 7 . Furthermore we make the following assumptions on 
the behaviour of d}i^ oia (x) / dx at endpoints of bands. Let z be a band end point. If the 
gap at z is a void, then 

1 da c ■ 

lim -^(x) = K, with < K < oo. 

x-+z ,xel ^j/\ x - Zq \ dx 

Similarly, if the gap at z is a saturated region, we suppose that 

]i m , 1 ( _ ^L( x ) \ K. W]1 h < K < -x. 



> z o, xei ^J\ x _ Zq \ \ c dx 
So the constraints are met like a square root. 

Remark. The one-point correlation function in the hexagonal tiling problem is shown to 
converge pointwise to the density of the corresponding equilibrium measure in [H Theorem 



3.12], see also Theorems 4.1, 4.2 and 4.3 In the (half-) hexagon tilings there is a single 
band corresponding to the intersection of a vertical line with the temperate zone, the 
surrounding gaps to the intersection of the line with the arctic regions. 

Finally we define the quantities involved in the asymptotic expressions for ~KN,k m a gap 
T and the band /, cf. [H Sect. 2.1.4]. The variational derivative of E c [fj] evaluated at 
/i = /i min is equal to 



"mm 



SE _ 5E c [pi,V] 

-{X) . 



8ji 5fi 



We have 



(x) = -2c / log \x - y\dfi c miQ (x) + cp(x). 
if x is in a void, 

if x is in a band, (3.10) 
if x is in a saturated region. 

The function ^f^ — l c defined in a gap V extends analytically from the interior. Furthermore, 

— r — / 

we denote by L c (z) and L c the function 




/ log|z-x|d/i min (x), 

J a 
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defined for z, in a gap T or a band I, respectively. L c (z) is analytic in z if $lz e T and Qz 
sufficiently small. L c has an analytic continuation to a neighbourhood of /. 



4 Statement of asymptotic results 

We have now introduced all the notation required to state our main results. The following 
results are verbatim copies of the corresponding results for the DOPE(2N,2k) ensemble 



in pQ, except for the "close to zero" assertions in Theorems 4.1| and 4.5 
The first results treat the band case and involve the sine kernel 



siri(7r(£ - v )) 
vr(£ - 77) 



(4.1) 



Extend S to the diagonal by setting S(£,£) := 1. Correlations for nodes close to are 
expressed in terms of the kernel given by 



sin(7r(£ - 77)) sm(vr(£ + rj)) 



Theorem 4.1 (Theorem 3.1 in [1]). Denote by Rm ,k \x\, . . . ,x m ) the m-point correlation 
function for the ensemble DOPE sym (N, k) with an even weight u>2N- Fix x 7^ in the 
interior of a band I := (a, (3) and let 



5{x) :-- 



c^$^(x) 



dx 



(4.2) 



Let . . . , belong to a fixed bounded set D cl in such a way that the points defined 
by 

jj)S(x) . 
xj :=x + Cn —, J = !,••• ,m 

are all nodes in X^. Consequently, Xj — >■ x as N — » 00. Then there is a constant Co,™ 
depending on D and m such that for sufficiently large N we have 



max 



5{x) 



tot (Sit®, $ 



l<i,jr'<m 



< 



c 



D,m 



N 



If (—(3,(3) is a band of dfi'^ ain (x)/dx, x = and the nodes are 

x r .= ti )6 ^,ti } >0,J = l,...,m, 
then there is a constant Cd,tu such that for sufficiently large N 

(1 \ m 
^— — J det^5 ,0 ^y,^ 



max 



l<ij'<m 



< 



Cd, 



N 
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Theorem 4.2 (Theorems 3.3 and 3.5 in [I]). For a fixed closed interval F in a gap T there 
are constants Kp and CV, m , such that for sufficiently large N and nodes X\, . . . , x m G FnY^ 

e -mK F N 

max \R^' k) (x!, . . . ,x m ) \ < C F , m 
xi,...,x m eF l\ m 

holds if T is a void. If T is a saturated region, we have the estimate 

-KpN 



max 

Xl,...,Xm&F 



\R^ k \x u ...,x m )-l\ <C F , 



N 



The result for nodes close to a band end point is expressed in terms of the Airy kernel 

Ai(£)Ai'fo) - Ai'(0Aifo) 



(4.3) 



Theorem 4.3 (Theorems 3.7 and 3.8 in [Tj). For each fixed M > 0, each integer m and a 
right band end point (3 separating I from a void, there is a constant G™(M) such that for 
sufficiently large N 



max 



X i 



detU(^,e. 



N , 



_(2AT)i/3 p o( /3 )_ 

where the max is taken over nodes xi, . . . ,x m G X N all satisfying 

f3 - M(2N)~ 2/3 < Xj </3- M(2N)- 1/2 , 



l<i,j<m 



< 



_/V(m+l)/3' 



the constant Bp is equal to (cf. Sect. 3.3) 



1 dn c 
Bp := lim -9^(x) > 



and £^ = {2NitcBp) 2 ^ (a;.,- — /?). is adjacent to a saturated region, then 



max 



(2iV) 1 /3 p o( / 5) 



N 



< 



iV 2 / 3 



where Hp l (M) is a constant, = (2Ntt(1 — c)BpY /A (xj — (3) and Bp is equal to 



2/3 



Bp = lim 



W y//3 - X 1 - C 



I p (x) _^nin (;E) 

c ax 



> 0. 



Remark. For DOPE(2N, 2k) this result is Theorem 3.7 and 3.8 in pp. Analogous results 
hold for left band end points. 

Another interesting statistic concerning particle systems is the fluctuation of extremal 
particles, e.g. the intersection point of a vertical line with the boundary of the arctic 
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region in random tilings. More generally, let B C Xn be a set of nodes and m an integer 
such that < m < min(#B, k). A well studied statistic is 

A^' k \B) : = P(there are exactly m particles in B) 
1 / d \ m 

= sj(-sJ det (I - ( K "*M 



t=i 



where is the operator on ^(Yn) with kernel K s ^(x,y) and K^k\ B its restriction to 

Let (3 be the rightmost band end point and (/3, b) be a void (resp. saturated region) 
Denote by x max the position of the rightmost particle (resp. hole, i.e. unoccupied node). 
Since the one-point function of DOPE(2N,2k) converges pointwise to cdfi'^ nin (x)/dx and 
DOPE sym (N, k) has the same asymptotic behaviour in regions bounded away from 0, 
one expects x max near the rightmost band edge (3. In this domain the correlation kernel 
approximates the Airy kernel. The latter kernel is the correlation kernel of the distribution 
of eigenvalues of a GUE matrix at the edge of the spectrum and the fluctuations of the 
largest eigenvalue are governed by the Tracy- Widom distribution v [T2] whose distribution 
function is equal to 

i/((oo,s]) = det (l-A\ [St00) ), 

where ^4|[ soo ) is the trace class operator on L 2 [s, oo) defined by the Airy kernel. The 
position rightmost particle (properly scaled) in DOPE(2N, 2k) is proved to be Tracy- 
Widom-distributed and the proof carries over verbatim to DOPE sym (N, k). 

Theorem 4.4 (Theorem 3.9 in [1]). Let the gap adjacent to (/3,b) be a void. Then the 
position of the rightmost particle x max in the ensemble DOPE sym (N, k) is described by the 
limiting distribution function 



lim P ( x max < (3 + S ) = det (l - A\ 



[s,oo) 



If the gap is a saturated region and x max describes the position of the rightmost unoccu- 
pied node (hole), then the same relation holds with c and Bp replaced by 1 — c and Bp, 
respectively. 

In DOPE sym (N, k) also the behaviour of the leftmost particle is of separate interest. Due 



to our assumptions in 3.3, lies in the interior of a band or a gap, the gap case being 



covered by Theorem 4.2 So we assume that no constraint is active in 0. We denote by S° 



the operator on ^(^) given by 



?0 _ sm Iww 7 ^ ~ sm [mm^ + .... 

7T(i-j) 7T(Z+J) 



where 8(x) is defined in equation (4.2). 
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Theorem 4.5. For the position of the leftmost particle x m \ n in the DOPE sym (N, k) en- 
semble we have the limiting distribution function 

^K" m, °-^) =det ( I_,s ° lB <»>)' 

where B(s) is the finite set of integers {0, 1, . . . , [s — 1/2J} and [s — 1/2J denotes the 
integer part of s — 1/2. 



5 Proofs of the asymptotic results 

5.1 Asymptotics for DOPE(N, k) 

We summarise weaker versions of some of the Theorems of pQ, which suffice for the esti- 
mates we need below. The first is on leading coefficients. 

Proposition 5.1 (Theorem 2.8 in pQ). The leading coefficients 7^ of PN,k( z ) an d jN,k-i 
of PN,k-i( z ) satisfy the asymptotic relations: 

J 2 N>k = C lNik e Nl ^(l + 0(l/N)) 

and 

ll, k - l =C^ k _ 1 e Nl ^(l + 0(l/N)) 
where C lNk _ 1 and C lNk are suitable constants and 7 remains bounded as N — > 00. 

Proposition 5.2 (Theorems 2.9, 2.13 and 2.15 in PQ). In suitable neighbourhoods of the 
different subintervals of [a,b] we have the following asymptotic representations of 7tN,k- 

1. Assume that J is a closed subinterval in a gap V . Then there is a neighbourhood Kj 
of J, and function A(z) analytic on Kj and uniformly bounded in N, such that 

knM = e^ (z)(2) (A(z) + 0(1/N)) 

holds. 

2. Assume that F is a closed subinterval in a band L Then there is a neighbourhood 
K F of F, and a sequence of analytic functions B N (z) defined on Kj and uniformly 
bounded in N, such that 

7r N , k (z) = e Ni;I ^ (B N (z) + 0(l/N)). 

3. Let z be a band endpoint. Then there is r > and sequence of functions Cn(z) 
analytic in \z — z \ < r and uniformly bounded for N — > 00, such that for \z — zq\ < r 
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More precisely we have the following results, the first of which treats the band case and 
involves the sine kernel (4.1). 

Lemma 5.3 (Lemma 7.13 in [I]). Fix x in the interior of a band I = (a, 0) and let 

-i 



6{x) :-- 



C : (X) 



dx 



Let £n,Vn belong to a fixed bounded set D cl in such a way that the points defined by 

. S(x) S(x) 
w := x + 4tv^^, z := x + r) N 



N 



N 



are both nodes in Xn- Consequently, w,z — >■ x as N — >■ oo. Then there is a constant Co 
depending on D such that for sufficiently large N we have 



max 



K N , k (w,z) - c d ^(x)S (Z N , m ) 



< 



C 



D 

w 



Furthermore, if Xn,i, %N,j — > x while i — j remains fixed, we have 

K N ,k (x N>i , x NJ ) = Sij(x) +0(l/N), 



(5.1) 



where 



sm 



Sij (x) 



C <Knin , 

p° (x) dx 



X) ■ TT(t-J i 



(5.2) 

□ 



Lemma 5.4 (Lemma 7.16 in [T]). Let (5 be the right band end point. For each fixed M > 
there is a constant Cp(M) > 0, such that for N sufficiently large 



max 



(ncB^ 



\2/3 



JVV3 



< 



Cp{M) 

AT2/3 



holds, where the max is taken over pairs of nodes x, y G Xn all satisfying 

f3 - MN~ 2/3 < ./-.// < I • MN- 1/2 , 



the constant Bp is equal to (cf. Sect. 3.3) 

and ( N = (NncB/s) 2 ^ (x - /9) and r] N = (NttcB^ 3 (y - 0). 
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Lemma 5.5 (cf. Lemma 7.4 in [T]). We have the exact formula for any node X — X jv n G 

X N 

1 



\/w N {x) 



7Tp°(x)N 



- h(x) e - |JV( ^ (x)-i c ) e - 1 (ATZ c +7) e -iV / a 6 log Ix-wld^Cy)^ § 



where the function 



Tjfiz) = cos 



'2ttN f*p°(x)dx y 



ttAT-1 / 



3 JV f^log\z—x\p°(x)dx 



is real analytic in (a, b) and bounded independently of N. 



Proof. Plug in the definitions and compare to the representation (3.2) of wn, this verifies 



the exponentials. The product in (3.2) is rewritten as 

N-l 



n 

m=0 



\%N,n ~ %N,r 



lim 

z— >XN 



Z — Xn , n 



COS 



2-kN p°(x)dx 



N-l/ \ 
m=0 \Z X N , m ) 



The limit of both fractions exists. De 1 'Hospital's rule applied to the first fraction yields 

1 1 



tc p°(x)N sin 
This completes the proof. 



2*Nf2fP(x)dx \ np°(x)N(-l) N - 



□ 



Remark. T^(z) has no poles at the nodes as the zeroes in the denominator are cancelled 
by the cosine. Furthermore the product in the denominator is asymptotically equal to the 
exponential for N — > oo. 

Proposition 5.6 (Lemma 7.12 and 7.14 in pQ). Denote by E the finite set 
{ao, (3o, . . . , ah, fih} of band end points. Let F be a fixed closed subset of the interval (a, b) 
such that F fl X N ^ and F fl E = ( and hence all band end points are bounded away 
from F). Then, for N sufficiently large, we have for all x,y G F fl X^ the estimate 

C 



(x — y)K Njk {x, y)e 2 >e 2 v s » w 



< 



H(N) 



(5.3) 



with H(N) = N and a constant C only depending on F. Furthermore, if x,y G X N are 
chosen from a sufficiently small neighbourhood G of E, the estimate holds with H(N) = 
TV 2 / 3 and the constant C only depending on G. 

Proof. Substitute the result of Lemma 5.5| and the assertions of Proposition 5.2 into the 
following formula (cf. (2.3)) 



(x - y)K Ntk (x,y) 

and estimate the uniformly bounded parts by a constant. □ 
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Remark. Recall that ^jjf-(z) — l c is strictly positive for z in the interior of a void and 
identically zero in a band. So, if at least one of x, y lies in a void, K^^(x, y) is exponentially 



small. Moreover, under the assumptions of Lemma 5.6, K^^{x^y) can be asymptotically 
non-zero only if x, y £ F D / and \x — y\ — 0(1/N) or, for x, y close to a band end point, 
\x-y\ = 0(1/N 2 / 3 ). 

We can now harness the results for DOPE(2N, 2k) for the ensemble DOPE sym (N, k). By 



equation (2.6) the correlation kernel for DOPE sym (N 1 k) is 

k n™( x i v) = K 2N,2k{x, y) - K 2N:2k (x, -y). 

We show that K 2 N,2k( x , ~y) is asymptotically neglegible and the results of the previous 
section also hold for DOPE sym (N, k). 

Lemma 5.7. Assume the situation for the DOPE sym (N,k) ensemble, i.e. [a,b] = [—b,b] 
and X 2N = —Yn U Y n for a set of nodes Y N = {x 2 N,N+j, J ' — 0, . . . , N — 1} and the weight 
function w 2 n is supposed to be even. This implies the same symmetry for the equilibrium 
measure /i^m- Moreover, we have the following estimates. 

1. Let x,y £ / n Y/v be a pair of nodes with x,y > 5\, where d~i is a strictly positive 
constant independent of N. There is a constant C depending only on 5i such that 

C 

\K 2N ,2k(x,-y)\ < — . (5.4) 

2. For nodes x, y £ J fl X 2 at, where J is a closed interval in a gap T, there are positive 
constants Di,D 2 depending on J, such that 

\K 2N , 2k (x,-y)\<D ie - D " N . 

3. Let M > and (3 be the right band end point. Then there is a constant Cm such that 
for any pair of nodes x, y with 

f3 - r(2N)- 2/3 <x,y </3 + r(2N)- 1/2 

we have 

C . 

(5.5) 

Proof. The symmetry of the weight w 2 n implies the symmetry of V 2 n and V (cf. equations 



(3.2) and (3.3)). Hence the external field y?(x) defined in (3.5) is even. It follows that 



the functional E c [/j] in (3.7) is invariant under the transformation fi(x) i— >■ /i(— x). By 



uniqueness the symmetry of /x^ in follows. 



The other assertions follow directly from Lemma 5.6, since in the first assertion \x — 



(—y)\ > 28\ is bounded away from zero. For the second assertion recall that by equation 



(]3.10|) ^f^{x) — l c is negative and bounded away from zero if x is bounded away from the 



band. Hence we have exponential decay in this case. For the third assertion notice that 
— (3 is the left band end point and \x — (— y)\ is asymptotically equal to 2/3. □ 
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Lemma 5.7 shows that, after an adjustment of the constants, the estimates for K2N,2k(x, y) 

"void" cases of Theorem 



in Lemmas 



5.3 



and 



5.4 



also hold for K s ^. Theorem 



4.1 



and the 



4.2 and Theorem now follow by simply taking determinants. For the "saturated re- 
gion" cases one has to take a detour over the dual ensemble, which describes the distri- 
bution of hole configurations, i.e. the probability p^ N ' N ~ k \y\, ■ ■ ■ ,yN-k) to find the sites 
yi, . . . , yN-k unoccupied. Saturated regions are voids for the dual ensemble and vice versa. 
The proofs are carried out for DOPE(N, k) in pQ and can be applied almost verbatim to 
DOPE s y m (N,k). 



6 Example: tilings of the half-hexagon 

The (Q, R, S)-hexagon is a hexagon with integer side lengths Q, R, S, Q, R, S and every 
angle equal to 120°. We study tilings thereof with 60° unit rhombi referred to as lozenges 
or simply tiles. The hexagon is filled without gaps and overlap, and no tile juts out beyond 
the boundary ("fixed boundary conditions"). The tiles occur in three different species 
(orientations) referred to as up-, vertical and down-tiles, see figure [TJ As we are interested 
in a certain symmetry class we restrict to hexagons with Q = 2k and R = S. To quantify 
things, we fix an ON-coordinate system and look at (symmetric) tilings of the (2k, R, R)- 
hexagon with corners (±V3R/2,±k) and (0, ±(k + R/2). The corners are numbered in 
counterclockwise order starting with Pi = (—\/3R/2, —k). If we focus upon the symmetry 
class of tilings w.r.t. the reflection in the x-axis, we can throw away the part below the 
x-axis and the chopped-in-half vertical tiles on the axis and obtain a tiling of the so-called 
(k, R) -half-hexagon, a model studied in [3]. 

Recall that tilings of the full (2k, R, i?)-hexagon map bijectively to families of 2k non- 
intersecting paths on the triangular point lattice 





L= \{— R --2) +q {--2) +r {---2)^- reZ 
with starting points in the set S = S + U S~ and end points in E = E + U E~, where 

5+ ={(-^'*+^)> < = 0...*-lJ, E+ = ^(^R, l + ^j, i = 0...k-l 

and S~ (resp. E~) denotes the reflection of S + (resp. E + ) in the x-axis. Admissible steps 
are (x,y) — > (x + \/3/2, y ± 1/2). To see this, connect in each up- and in each down-tile 
the midpoints of the vertical sides by a straight line segment (the decoration depicted in 
figure [T]) . 

In the same fashion tilings of the half-hexagon are mapped to families of k non- 
intersecting paths with starting points in S + and end points in E + which do not touch the 
x-axis. Denote by L m , m — 0, . . . , 2R the set of numbers such that {(—R + m) a/3/2} x L m 
is the intersection of the vertical line x = (-R + m)\/3/2 with the lattice L and the 
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Figure 1: Tiling of a (2k, R, i?)-hexagon (without symmetry), tiling of a (k, i2)-half-hexagon, 
up-, vertical and down-tiles 



(2k, R, i?)-hexagon. Note that by symmetry L m = L 2 ^ m and, for m < R, \L m \ = 2k + m. 
Then {(-R + m) \/3/2} the set of possible points where a family of lattice paths 

corresponding to a tiling can intersect a vertical line after m steps. Denote by the set 
of positive elements of L m . 

Proposition 6.1 (Theorem 4.1 in [S], Lemma 2.2 in |4j). Consider the sets of families 
of 2k non-intersecting lattice paths with starting points in S and end points in E (tilings 
of the (2k, R, R) -hexagon) and of families of k such lattice paths with starting points in 
S + and end points in E + not touching the x-axis ((k, R) -half -hexagon) to be equipped with 
their respective uniform distributions. 

1. Let X\ < x 2 < . . . < x 2 k, Xi £ L m for i — 1, . . . , 2p. Then the probability of a family of 
lattice paths to intersect the vertical line x = (-R + m)\/3/2 at ordinates X\, . . . ,x 2 k 
is equal to 

1 2p 

P m (xi,x 2 , . . . ,x 2k ) = ~~- \\w(xi) Yl (xj-xi) 2 . (6.1) 

i=1 l<i<i<2p 

2. Let Xi < x 2 < . . . < xjz, Xi G for i = 1, . . . , k. Then the probability of a family of p 
lattice paths not touching the x-axis to intersect the vertical line x = (-R + m)\f?>/2 
at ordinates X\, . . . , x^ is equal to 

c^^",^)=in*) n w-*?) 2 - ( 6 - 2 ) 

m i=l l<i<j<P 
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In both cases the weight function w = w m is even and it is equal to 

W ^ = (m/2 + k - l/2± z)\(R - m/2 + k - 1/2 ± z)V ^ 
Z m and Z^ m are normalisation constants. 

Note that every point in {(— R + m)\f?>/2) x L m where none of the paths intersect is 



a centre point of a vertical tile. So the dual ensemble of (6.1) describes the distribution of 
vertical tiles along the line x = (-R + m)\/3/2. 

The weight function w turns out to belong to the Associated Hahn family. 
6.1 The Hahn and Associated Hahn ensemble 

Let [a, b] = [-1/2, 1/2] and X N = {x N>n = -1/2 + (2n + 1)/2N, n = 0, . . . , N - 1}. The 
discrete orthogonal polynomial ensemble with weight function 

.AHE ' 



(x N , n ; P, Q) 



n\{P - 1 + n)\(N - 1 + n)\(Q - 1 + JV - 1 - n)\ 



is called Associated Hahn Ensemble (AHE) [H |6] with parameters P and Q. Its dual 
ensemble is known as Hahn Ensemble (HE) with weight w^ E (-;P, Q). In [lj Sect. 2.4.2] 
the equilibrium measure /i^ in for the family w% E ( ■ ; AN + 1, BN + 1), A, B > 0, c 6 (0, 1) 
fixed is computed. It turns out that there is exactly one band interval. For A = B it is an 
interval (— /3, (5) enclosed by two gaps of the same type. If c < ca = V 'A 2 + A — A those 
two gaps are voids for HE (and hence saturated regions for AHE with c substituted by 
1 — c) and if c > they are saturated regions for HE (voids for AHE, c substituted by 
1 — c). For the right band end point one has 

a = y / c(l-c)(2A + c)(2A + c+l) 

P 2(A + c) ' { ' ' 

The case c = ca is exceptional, leading to gaps consisting of single points. 

6.2 The parameters in the (2k, R, i?)-hexagon 

With the change of variables z = n — (m/2 + k — 1/2) and the substitutions 2k + m = N 
and R — m + l = Pin the weight function w(z) in ( |6.3 ), we see that w is from the AHE 



family with the parameters P and Q both equal to R — m + 1. With view on the half- 
hexagon and the general results of the previous sections we consider the parameter m in 



the distributions in Prop. 6.1 to be even, m = 2p, since this implies that the set L 2p has 
an even number 2N = 2k + 2p of elements. For m — 2p + 1, we have G L m and the 
remarks at the end of Section |2] allow us to treat that case in the same fashion as below. 

We choose R such that R/k — > A as k — > oo. We scale the (2k, R, i?)-hexagon by 1/k, 
such that in the limit its intersection with he x-axis is the interval [-y/3\/2, V3A/2]. We 
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want to study the intersection behaviour of the rescaled family of 2k lattice paths with a 
line x = r E (— \/3A/2, 0), i.e. the parameter m = 2p shall satisfy the asymptotic relation 
— >■ t for — >■ oo and hence 

2p/A; -»■ 2t/V3 + X. 

This implies for the ratio c of the numbers of paths ( "particles" ) 2/c and nodes 2 N 

2k/2N = 2k/(2k + 2p) ->■ c := = . 

7 /V 7 2 + A + 2r/V3\/2 

Moreover for the parameters P = R — 2p + 1 and A we have 

(P - I) /2k = 2NA/2k -> A + r/v 7 ^ 

and hence we can choose 

A = (-r/V3)c. 

The set of nodes X 2 n = {l/k, I G -^2p} lies for large k in the interior of the intervall 
[— 1/c, 1/c]. We can compute the band endpoints according to (6.4). Bearing in mind that 
that we are dealing with the dual ensemble of the Hahn ensemble, we have to substitute c 
by 1 — c. After a proper rescaling by 2/c the band endpoints read 



±-P = ±V\ + MI 1 



Vv^A 



T 



r e 



v^A v^A 



(6.5) 



The points with these ordinates and abscissa r lie precisely on the boundary of the inscribed 
ellipse of the rescaled hexagon. Th e neighbou ring gap is a saturated region if r < r := 



2(2+0' 



i.e. 



c > c A , 



cf. 



6.1 



Theorem 



4.2 



yields that above the inscribed ellipse the line x 



for the one-point correlation function 
= t is fully packed with intersection 



points (and hence with the blue "up-tiles"), with probability exponentially close to one. If 
< t < 0, i.e. 1 — c < ca, then the gaps are voids and the part of the line x = r 



2(2+0 



above the ellipse is devoid of intersection points and hence fully packed with vertical (red) 
tiles with probability exponentially close to one. Close to a point with ordinate > in 



the interior of the ellipse, by Theorem AA_ the one-point correlation function takes a value 
strictly between and 1 and hence both intersection points and vertical tiles occur with 
positive probability. This reflects the "arctic phenomenon" present in the large tiling of the 
half- hexagon depicted in Figure [2j namely the tiling being highly ordered at the corners and 
unordered in the interior, with a sharp transition along the inscribed ellipse. Furthermore, 
the fluctuations of the highest vertical tile (r < tq) resp. the highest intersection point 
(r > Tq) are governed by the Tracy- Widom distribution by Theorem 4.3 Close to the 



x-axis approximations of the correlation functions for the intersection points of the paths 
with the line x = r are also given by Theorem AA_ in terms of the kernel S°. This is in good 
accordance with a result in [13J on non-intersecting Brownian excursions, the continuous 
counterparts of our path model corresponding to tilings of the half-hexagon. 
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Figure 2: Arctic phenomenon in a (32, 64)-half- hexagon 



7 Conclusion 

We have applied the results of the monograph pQ to a general discrete ensemble a special 
case of which, namely random tilings of the half-hexagon, is applied in the study of minors 
of certain random matrices [I]. Our results on the correlations of tiles on a vertical line 
and the "arctic-half-ellipse" result complement the work in on the half-hexagon model. 

The continuous counterpart of our general ensemble appears in the study of non- 
intersecting Brownian excursions [T3] which matches the lattice paths formulation of the 
half-hexagon model. In [13] also the distribution area below the lowest path is addressed 
and the first moment thereof is computed. From a combinatorial point of view it would 
be interesting to (asymptotically) compute the expected area below the lowest path in a 
family of lattice paths. Furthermore, higher moments or the limiting distribution of area 
would be interesting if, say the number of paths is fixed and their length tends to infinity. 
This has so far only been fully accomplished for a single path in l l \]. 
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